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1 Two-dimensional contour plots showing regions within which, f^ , the maxi- 
mum electron gyrofrequency at the mirror field B^, exceeds 27, 20, 16 or 
10 MHz (shaded areas and dashed line). The contours are projected onto the 

= altitude of the cloud tops. Also drawn are isointensity contours in 1G incre- 
ments oi'Tupiter’s surface field. Footprints on the surface are plotted for 
L = 6 (^”), L = 5 (“X”), L = 4 (“*”), L = 3 (“0”) and L = 2 (“#”) field lines. 

The azimuthal spacing of the field lines in the equatorial plane is 10°. Con- 
jugate points at any longitude can be found by noting that in the two hemi- 
spheres, near Xjjj = 0°, footprints denoted “S’* are conjugatefc ; (04 model) . . 26 

2 Similar to Figure 1 , but using the JPL model 27 

3 A plot of i versus Am at 27 MHz, using the 0 4 model with = 19°. Gaps 
in the curves for L = 3 and 5 result from localized regions where f M is less 
than 27 MHz. Curves on the right are solutions of Aj.jj ~ Xjjj' + n, while 
those to the left are obtained from Ajjj = Xjjj' - fi. The curves represent 
the central meridian longitude directions of the intersection lines of the 
emission cone with the ecliptic. Note that there is no longitude for L = 2 at 
which f M > 27 MHz (cf. Figure 1). Preferential beaming into the ecliptic — 
small values of r — occurs at the positions of the Io independent Early and 

Main Sources (B and A, respectively) 28 
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' 4 Similar to Figure 3, but for 20 MHz. For clarity, curves are not plotted for 
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(Figures 4a and 4b, respectively). In the north, minima in i (peaks on the 
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(C) longitudes. See the text for a discussion of the southern hemisphere 
sources . . .... 29 

5 Similar to Figure 4, but using the JPL model. Note that for L = 2, the 
allowed range of Am is much less than it was using the 0 4 model (cf. Fig- 
ure 3). In addition, at L = 6 in the north (Figure 5a), small values of i do 
not extend to non-Io C longitudes. See the text for a discussion of the 
southern hemisphere source locations (Figure 5b). . . . . . . . , . . , 30 

i. A plot of t versus Am at 10 MHz, using the 0 4 model. Use of the JPL 

,/ model produces identical curves. In the north (Figure 6a) \p = 40°- in the 
south (Figure 6b) 'P = 60°. See the text for a discussion of the source 
locations 31 
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ABSTRACT 

A geometrical model is presented in which the apparent source 
locations of the lo-independent decameter radiation are computed. 
The calculations assume that the radiation is produced by stably 
trapped electrons radiating near the local electron gyrofrequency 
and that the emission is then beamed onto a conical surface. The 
maximum occurrence probability of noise storms is associated 
with regions in the Jovian magnetosphere where the axis of the 
emission cone is most inclined toward the Jovian equatorial plane. 
The calculations utilize and compare two of the octupole spherical 
harmonic expansions of the Jovian magnetic field constructed 
from data accumulated by the fluxgatc and vector helium magne- 
tometers on board Pioneer 1 1 . 


A BEAMING MODEL OF THE IO-INDEPENDENT JOVIAN 
DECAMETER RADIATION BASED ON MULTIPOLE MODELS 
OF THE JOVIAN MAGNETIC FIELD 
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I. INTRODUCTION 

Jupiter’s low frequency radio emission at decameter (DAM) wavelengths has two 
components, In the first, the occurrence probability of a noise storm depends on both 
the Jovian central meridian longitude (CML) that faces Earth and the orbital phase of 
the innermost Galilean satellite, lo. Secondly there is a significant probability of observ- 
ing DAM radiation at particular central meridian longitudes, independent of Io control. 
(For a recent detailed review of the observations the reader is referred to Carr and Desch 
(1976).] In this paper, we concentrate on this Io-independent component of the decam- 
eter radiation, adopting the viewpoint that although the origins of the Io-controIIed and 
Io-independent radiation are both profoundly influenced by the local properties of the 
Jovian magnetic field, the physical emission mechanisms and source locations of the two 
components can be distinct. 

The Jovian magnetic field determines the source locations of the radiation, and 
the magnetic field measurements made in the inner magnetosphere on Pioneer 1 1 have 
greatly increased our knowledge of the detailed structure of that field. Models have been 
constructed which permif extrapolations of the field down to the Jovian cloud tops. The 
measurements indicate that within six Jovian radii (Ry), models of the Jovian magnetic 
field will reasonably fit the data if they include internal source terms up to the octupole. 
The surface fields that one infers in the northern hemisphere reach 14G, while in the 
south the maximum is approximately 10G. In a 14G field, the local electron gyrofre- 
quency is 39 Mliz — very close to the maximum observed frequency of the Io-controlled 

f :V 
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DAM component (Warwick 1964). The highest electron gyrofrequency in the southern 
hemisphere is 29 MHz — close to the maximum of the lo-independent component (Wilson, 
Warwick, and Libby 1968; Thieman and Smith 1978). We shall assume that this 29 MHz 
cutoff in the lo-independent component is a physical effect not an observational limita- 
tion treated by the combined effects of limited antenna system sensitivity and the greater 
. \ 

intensity of the Io-controlled component, as was suggested recently by Desch (1976). 

\V' 

Based on this evidence, we conclude that the decameter radiation is produced near 
the local electron gyrofrequency and, furthermore, that emission associated with the lo- 

/r . 

independent component originates from electrons stably trapped in the Jovian magneto- 

sphere, for trapped particles cannot reach fields higher than the maximum in the south. 

Goldstein and Eviatar (1972) proposed that trapped electrons could produce the Io-inde- 

im- 
pendent DAM via a cyclotron instability driven by a loss-cone distribution. However in 

that paper, a centered dipole field was assumed with a field strength far in excess of what 


we now know to exist. Thus, clearly our increased knowledge demands some modifica-_. 


tion of those ideas. Using the high-order multipole models of the magnetic field to con- 
struct the beaming model developed below, we can now estimate the longitudinal locations 
of che radiation sources in the Jovian magnetosphere. From these source positions, the 

longitudes from which radiation is beamed toward Earth can be found, and one canrpro- 
& ' j? 

ceed to develop more complete theories of the emission mechanism. Suih a theory is 


proposed in the accompanying paper (Goldstein and Eviatar 1978). 




Implicit in our analysis is the assumption that the DAM escapes in the R-X mode, 


\\ 


near the local electron cyclotron frequency. In fact, we will be most interested in the 


free escape propagation band (Band 111) for which the angular frequency co exceeds cor, 


( - > 
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where (Stix 1 962) 

cj r « wln e i[i + (» + 4w pe 2 /n e 2 ) w l ^ (l) 

in which co pe is the local electron plasma frequency, and £l e = -|e|B/Mc is the electron 
Larmor frequency. The emission frequency will then be close to the gyrofrequency as 
lon/as w pe « |S2 e l- This implies a local electron density, N e , less than 10 6 cm 3 , which 
seems likely everywhere in and above the Jovian ionosphere (Fjeldbo, et al. 1976; and 
Atreya and Donahue 1976). 

In the development presented below we examine the implications of the assump- 
tion that stably trapped electrons are the exciters of the Io-independent DAM. This can 
be done regardless of the particular instability mechanism. 


II. TRAPPED ELECTRONS IN A MULTI POLE FIELD 
a) Source Locations 

In our analysis we have utilized several models of the Jovian magnetic field to 
examine the sensitivity of our conclusions to the different extrapolations of Pioneer 1 1 
data to the very low altitudes where 10-30 MHz gyrofrequencics are encountered. The 
models we have used include the O4 model of Acuna and Ness (1976), containing tiftce 
internal (1) and no external (E) sources, and the PI 1 (I3-E0) and PI 1(I3-E1) models pro- 
posed by Smith, Davis, and Jones (1976). In the last, because we wished to utilize 
a time-independent model in a frame rotating with the interior sources, we discarded the 
g] * andTTj * exterior source coefficients as recommended by Smith et ah, leaving only 
the "gj 0 exterior coefficients. Surface maps of the magnetic field computed from all the 
<r models are very similar, especially in the northern hemisphere. The most significant dif- 
ference is the strong tendency of the Smith et ah models to place' the highest southern 
hemisphere fields at higher System III (1965) longitude, Xm, than does the O4 model 
(330° compared t& 255° in the 0 4 model). As we shall see, the accessibility of trapped 

V 

electrons to regions 6f high gyrofrequency is quite sensitive to the relative locations of 
the high field regions in the t%vo hemispheres. Consequently, we have found it necessary 
to present some of our conclusions in terms of both the 0 4 and PI 1 (13-El) models, (In 
general, the PI l(I3-E0) model is very similar in most respects to PI 1(I3-E1)). To simplify 
further reference to the various Smith et ah models, we will denote the PI 1 (13-El ) model 
as the “JPL” model, and continue to use “O4” to designate the Acuna and Ness model. 

The locations in the Jovian magnetosphere that give rise to a particular gyrofre- 
quency are found by integrating the equations describing the trace of a magnetic field 






line in spherical coordinates. The integration begins at some point (r,0) in the zenographic 
equatorial plane and proceeds to calculate B(r,0,tf>) down to the cloud tops in the northern 
and southern hemispheres. The maximum mirror field, B^j , associated with a particular 
field Jine is then the lower of the magnetic field values at one of the two ends of the 
field line. [In locating the foot of a field line, one must include the planetary oblateness. 
For ♦he radial distance of the clout . ?s, we have used (Kahle, Kern, and Vestine 1964) 
} R/R„ = (1 + ecos2 0 )-K (2) 

with e = 0.14371 (Anderson, Null and Wong 1974), where 0 is colatitude and R^ = 
71,372 km.] The gyrofrequency, fj^j, corresponding to By, is then the highest that 
can be sampled by a trapped electron on that field line. 

If one evaluates f M , beginning the field line integration at various equatorial 
radii and azimuths, one can delineate the regions within which radiation above some 
predetermined frequency originates. Given a specific field njodel (e.g, O 4 or JPL), the 
possible source locations of the decameter radiation can be found as a function of fre- 
quency: higher frequencies coming from increasingly restricted areas and lower alti- ; 

'■/ < 

tudes within the magnetosphere. In our calculations, we rsive found it sufficient to 
consider equatorial radii L < 6 . Here, L denotes a magnetic field line that intersects 
the equatorial plane at r/R u s L. The analysis can be readily extended to larger L- 
values, if desired. [The reader will note that L has been defined with respect to the 
zenographic, rather than the zenomagnetic, equator for computational convenience. 

Our analysis is insensitive to small changes in L that would; result if we used the more 
customary definition,] ./ - 


Figure 1 illustrates the results of such a calculation using the 0 4 model. Shown 
in a two-dimensional polar plot are the projections onto the surface of regions within 
which f\j exceeds 10, 16, 20 and 27 MHz in the northern hemisphere (Figure la), and 
southern hemisphere (Figure lb). Only those regions accessible to electrons on L-shells 
of 2, 3, 4, 5, and 6 are indicated. Also plotted in the figure are the footprints at the 
cloud tops of those same L-shells. Finally, the figure includes contours of constant 
surface magnetic field intensity plotted in one-Gauss increments. The maximum gyro- 
frequency sampled by trapped electrons at altitudes above the Jovian cloud tops is 
found to be nearly 29 MHz — close to the maximum of the lo-independent DAM. 

Several features on the {.jots are noteworthy. First, it is clear that f M exceeds 
27 MHz over a small region, very confined in azimuth, accessible to electrons on L- 
shelh between three and six. In the northern hemisphere, the f^ > 27 MHz region 
occurs at lower latitudes than the corresponding region in. ,the southern hemisphere. As 
the mirror point gyrofrequency decreases, larger emission regions become accessible. By 
16 MHz, .particles on the L = 6 flux tube can radiate from all azimuths. As the mirror 

| j 

point gyrlifrequency reaches 10 MHz, electrons on L-shells > 2 can radiate at all azi- 

!! 

mut^s (except for a small region near Xjjj = 280°). This is not to say, however, that 
radiation at these low frequencies should be observable from all CML. The detectabil- 
ity of the radiation depends critically, of course, on the excitation mechanism, beaming 

/{. 'J 

pattern, and on propagation effects such as refraction. Also, it is not known a priori 
which L-shells are populated with electrons that produce DAM. Although it is gen- 
erally assumed that the lo-contrqlled component originates at L = 6 (the orbit of lo), 
no such general statement can be made about the component independent of Io control. 
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For comparison, in Figure 2, a similar plot has been constructed using the JPL 
model. As mentioned above, the most significant difference in this model is the location 
of the maximum southern hemisphere field at Am = 330°. As a consequence, the loca- 
tion of regions within which trapped electrons sample gyrofrequencies in excess of 27 
MHz is now a very small area between L — 5-6, at higher Ajj| than ill the O4 model. An- 
other difference between the models is that the highest surface field in the JPL model is 
greater than 15G, but is about 14G in the O4 model. [It should be noted that this model 
seems unique in this respect. The other Smith et ah (1976) models we have examined, 
as well as that of Smith ct ah (1975), predict maximum fields of about I4G, but still 
systematically higher than those found in either the O4 model or the unpublished Acuna 
and Ness model O5, which contains the same number of coefficients as the JPL model 
that we are using.] It is also noteworthy that the highest gyrofrequency that can be 
sampled at all azimuths is 16 MHz in the 0 4 model and 20 MHz in the JPL model. This 
difference would also hold true for precipitating electrons. 

b) Radiation Pattern 

/if, as we suppose, the decameter radiation is generated by an instability of an elec- 
■/ • 
tron in gyroresonance with a wave having to near co R , then the growth rate of such a 

wave will have a maximum in a direction determined by both plasma kinetic processes 
and by large-scale gradients in the local magnetic field. Regardless of the instability mech- 
anism, there is evidence that the decameter radiation above 10 MHz is beamed into a 
conical surface with a fairly wide half-angle, 'P (v. Dulk 1967, and Carr and Desch 1976) 
which is about 80° if L = 6 is the flux tube of the source. 
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Using either 'Magnetic Reid models, it is possible to find the apparent source loca- 
tions in CML from which radiation would be received in the ecliptic from the !/#>;• sides 
of the emission cone. In general, the axis of the cone will be inclined to the equator by 
an angle t, while the cone will intersect the equator along a curve whose asymptotes are 
two straight lines intersecting at an angle in. Once fi is determined, tire source locations 
are easily computed. From Goldreich and Lynden-Bell (196*)). p can be found from 

tan fj = ^/(taii- 'l' cos- t - sin- R 


(3) 


'flic angle t is defined by t = cos -1 (B c /B), where B c is the projection of B onto the equa- 
tor. If we describe the Jovian field temporarily in terms of a conventional right-handed 
spherical coordinate system in which 0 is Hast Longitude, then B c is given by 

B c — \/[B^- + (B r sin Q + cost?)-] 


.,) u c ~ V l u ti~ ' V u r iul1 u ' W9UJ-I (4) 

Defining g » sin -1 (B^/B c ), the System 111 longitude toward which B c (-B c ) points in the 
northern (southern) hemisphere is X m * = <? m -p\ (Because System 111 is a West: Longitude 
convention, 0 jjj = In - $). Finally, the appoint source locations seen from the equator, 
arc Am(±p) = Xu/ ± P> with plus and minus i/,bsig?,a ting the two sides of the emission 
cone. The apparent source locations are thus seen to be rather sensitive functions of 4f 
_ and i. 

Once the emission cone angle, 4' , is chosen, the apparent source locations can be 
found for emission generated at any point on a magnetic flux tube. When i is large, the 
emission cone may not intersect the ecliptic at all, while for small values of t, the cone 


|K>ints nearly directly into the ecliptic, resulting in source locutions separated by 2(i a 24 f . 


‘-Whether or not niditjion at some frequency, f, will actually be observed in the 
ecliptic from a particular CML will depend on severat “factors. First, of course, f cannot 






exceed f M at the apex of the emission cone; and second, given 'I', the emission cone must \ 
intersect the ecliptic, he. i < But even these conditions may not suffice. The occur- 
rence probability, P, of the various decameter sources appears to be a very sensitive func- 
tion of i. This has been inferred indirectly from the large variation of P with Dg. [Dg 

p 

is the Joviceniric declination of the earths the angle of the Jupiter-Earth line with re- 

0 

spect to the plane of Jupiter’s equator, varying by ± 3.3° during a 1 1.9 year Jovian orbital 
period. We will return to a discussion of the Declination Effect below.J The evidence 
suggests that P is largest when the magnetic field in the emitting region is most inclined 


h , 

toward the earth (v., ej*. Dulk 1967, and Schatten and Ness 1971). We fire thus led to 
the following assumption: along a given L-shell, the highest occurrence probabilities are 

■> •- s " V; 

located at minima in plow of i versus Ajij. Furthermore, there should be a qualitative 


resemblance between i plotted against Am, and plots of P versus CML. IV 

A quantitative justification of this assumption must await a more detailed treat- 
ment (including refraction effects) of the emission mechanism. We now proceed to con- 
sider in some detail the variation in i with CML, 'P, L and hemisphere, and its relationship 
to observed DAM morphology. 


c) Location and Occurrence Probability of the lo-Independent DAM 
To find the regions of CML where radiation is most strongly beamed into the 

;/ 

ecliptic (minima in a plot of i versus Am), it is necessary to choose 'P. Dulk (1967), in 

j! >rv 

analyzing the lo-related components, concluded that 'P = 79° woiild associate the Io-A 
and B (Main and Early) sources with two sides of a single emission cone. Goldreich and 
Lynden-Bell (1969) found support for this value from independent theoretical arguments. 
However, as there is no a priori reason for the non-Io-components to be generated by the 



same mechanism, we treat ^ for the moment as a free parameter. The angle will be cho- 
sen to give the observed separation between the non -Io A and B sources in the northern 
hemisphere. At 10 MHz we have found it necessary to choose a value for 'I' in the south- 
ern hemisphere which differs from the one that gilies the best fit for the northern hemi- 


sphere Sources. ^ 

Desch et ah (1975) and Desch (1976) have discussed the distinction between the 
lo-controlled and Io-independent sources at various frequencies. They noted that at any 
given frequency, the Io-independent radiation is less intense than the Io-controlled. Thus, 
for example, at 26.3 MHz, DAM noise storms having peak fluxes of less than 2 x 10 4 Jy 
(1 Jy = 10 ~ 26 W m ~ 2 Hz -1 ) are predominantly independent of Io control. At that fre- 
quency, high occurrence probabilities are centered near 1 10° and 250° CML; non-Io 
B and A, respectively. =\-. 

-'y ' 

We begin oUt interpretive discussions at high frequencies (27 MHz) where the 
analysis is relatively simple, and then progress by stages to the more complicated situation 

at 10 MHz. A gyrofrequency of 27 MHz (and by implication, the apex of the emission 

({. () 

cone) is found over only a very limited range of X|jj (cf. Figures 1 and 2). Thus the two 

- ■ // ■ '7 

sides of the emission cone will only intercept the ecliptic over a limited range of Am* 

(1 

^^hen the C 4 model is used, the longitude range witliin which one can obtain f = 
27 MHz in the' {northern hemisphere coincides with the greatest inclination of the emission 
cone toward the equator (he., 1 attains its minimum values). In the south the conjugate 
longitude range occurs near the maximum values of t. Hence, even before choosing a 


} 

value for 'P, it is clear that emission at 27 MHz will be most easily beamed into the equa- 

■ ■■■ '' V 

tor from the northern hemisphere. This is encouraging because, as mentioned above, Vv 
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studies of the Io-controllcti sources have led to the conclusion that the source locations 
are correlated with small values of t. 

In choosing it was natural to try first the value of 19° obtained by Dulk (1967). 
The results are shown in Figure 3, where t is plotted against Am for various L-values and 
for both hemispheres. Note that there is no curve for L = 2 because f M < 27 MHz every- 
where on that L-shell. The two symmetric curves (Am = X.ju ± p) stand out clearly. 
Because i in the southern hemisphere is large, radiation from there would not easily reach 
the ecliptic; nor is there observational evidence fti; left-polarized radiation at 27 MHz. 
(Electromagnetic waves propagating in the R-X mode would be received at Earth as right- 


p Jlarized if emitted from the northern hemisphere and left-polarized if emitted from the" 


south.) In the northern hemisphere, with Sk = 79 u , the two rides of the emission cone 

along L = 6 intercept the equator at Am centered near 120° and 270°. This corresponds 

\ 

closely to the observed locations of non-Io B and A, respectively. When L-values less than o 
six are considered 'I' as 60° fits somewhat better. In general, it will not be possible to 
use this beaming model to determine the L-shell on which lo-independent radiation 
originates. 

Thus, with 'I' = 19° (or 60° for L < 6), minima in i are located where Ajjj equals 
the observed CML of the non-Io A and B sources. However, as expected, the resemblance 
of Figure 3 to observations of occurrence probability is only qualitative. In particular, 

V , '' "./■ 

the observed occurrence probabilities of Sources A and B are never equal, and generally / 
non-Io A occurs more frequently than non-Io B. Whether this asymmetry may be a con- 
sequence of refraction is a question to which we return below. 



A^imilar plot can be constructed using the JPL model. However, at 27 MHz, only 
L values otV5 and 6 allow 5* 27 MHz, and then only over a very narrow range of lon- 


gitudes (cf. Figiiie 2). In fact, i ^ 50° in both hemispheres. Thus, in this modei_and at 


A. 




this frequency, radiation is not readily ^beamed into the ecliptic at any longitude. 

In general, as f M is decreased and as the altitude of the source region increases, 
the differences between the O4 and JPL models become smaller, disappearing altogether 
by 1 0 MHz. These discrepancies probably result from the large extrapolation of the Pio- 
neer 1 1 data from the region where they were obtained to the very low altitudes at which 
fj^j ~ 20 MHz. Gyrofrlq uencies of 27 MHz are found only close to the cloud tops, and 


the octupoie contributions, which are the least well-determined, are very important, while 
lower gyrofrequencies (~I0 MHz), are found at considerably higher altitudes where the 
extrapolations of the Pioneer 1 1 data are not as severe. Thus it is not surprising that the 

jf 

magnetic field models tend to become similar at lower frequencies. One should also keep 

to 

in mind that 27 MHz is very close to the observed upper limit of the non-feBAM. Io- 
independent radiation is clearly present at ,26.3 MHz (Desch^et ah 1975), and has proba- 
bly been seen at 27.6 MHz (Thieman and Smith 1978), but thus far no definitive identi- 
fications have been made at higher frequencies (v. Wilson, Warwick, and Libby 1968). 

At 20 MHz, - 79° again appears to locate the minima in 1 versus Ajjj closer to 
the observed source locations of non-Io A and B, than does, say ip = 60°, at least for 

a ' 

radiation originating near L ££ 5-6. In Figure 4, the computed values of 1 are shown using 
the O4 model. For clarity, northern and southern hemispheres are plotted separately. In 
the northern hemisphere, broad minima in t occur near Aju = 80 °- 150 ° and 240 °— 330 ° 
at L as 6. Lower values of 1 occur at lower L-values, again suggesting the possibility Jhiit 

V\ 
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the non-Io DAM may originate inside a large region of the magnetosphere bounded by L 
between 2 and 6 . If the radiation is dominated by sources at low L-vulues, then as 60° 
provides a better fit for the location of non-Io A and B. There is no way to distinguish 
between these, two possibilities Within the context of this model. 

In the southern hemisphere, narrow minima in i occur near Aju as 15° and 160°. 
The narrowness of these minima suggests that it might be comparatively difficult for rad- 
iation excited in the southern hemisphere to reach the ecliptic plane. From L - 6 there 
is an additional area of small i near Am =* 90° and 300°. 

The works of Dcsch (1976), Bozyan and Douglas (1976), and Thieman (1977) 
contain extensive analyses of DAM observations at 18 or 20 MHz. The non-Io controlled 
components are most likely to be observed near 160° CML (non-lo B), 250° CML (non- 
lo A) and 32Q° CML (non-Io C or Late Source). As before, the Early and Main Sources 
(B and A) have a natural explanation in terms of this beaming model. In addition, based 
on the O 4 model, the radiation from Source C may also originate in the northern hemi- 
sphere, because small values of 1 extend to 330° on the L = 6 plot. However, Kennedy 
(1969) has reported that Source C is often left-circularly polarized, which suggests that it 

■r 

originates in the southern hemisphere. Tire southern hemisphere minimum on L = 6 at 
300° in Figure 4b is close to Source C, while the second minimum near 90° could be an 
indication of an Io-independent Source D, 

Support for the idea that Source C comes from the southern hemisphere can be 
found from the JPL model (Figure S). The northern hemisphere minima in 1 (Figure 5a) 
are similar, though narrower than those in Figure 4a. In particular, for the region near 
Ajii — 240°, the minimum in t no longer extends into the Source C region. [Also note 
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that < 20 MH? virtually everywhere on L - 2.] In this model, the southern hemi - 1 
sphere t minima are found near Am a i30° and 330°. Emission beamed lrofii 130° CML 
in the southern hemisphere could be masked by emission from the same CML in the 

north. Emission beamed from the minimum in t at Am — 330° in Figure 5b is within 

\y 

the observed non-Io C source. It is also interesting to note that non-Io C is not observed ° 
above 22 MHz; and in neither of the magnetic field models is there any possibility of 
beaming into the ecliptic from that source region in either hemisphere at 27 MHz. In fact, 
there is a general tendency in the southern hemisphere for the highest frequency source 


regions (f^j > 20 MHz) to be located where * is largo (cf. Figure 3). fn contrast, in the 


north the highest frequency source regions are located where i is small. 

The two multipole magnetic field models both suggest an enhanced probability 
that radiation emitted into a 79° (L s 6) or 60° (L as 3) emission cone angle will be seen 
in the ecliptic as coming from the observed source locations of non-lc A and B. Concern- 
ing Source C, the information provided by these models is ambiguous. Taken together, 
they argue for non-lo C coming from the southern hemisphere, with some contribution 
irom an extended Source A in the northern hemisphere. However, while this geometrical 
construction provides a clear framework within which the northern hemisphere sources can 
be understood, the location of Source C is not as well determined. The reasons for this 
asymmetry are considered in detail in the accompanying paper by Goldstein and Eviatar 
(1978), but we will mention one difficulty here. Because there are so many similarities 
between the to-controlled and lo-independent DAM sources, it is natural to assume that 


regardless of the specific excitation niechanisnjfihe source locations should be closely re- 


f- 


lated. The Io-C source is observed when The phase of lo is near 240° as measured from 

if 


superior geocentric conjunction. This places the lo-C source on the +fi side of the emis- 

!/ 

sion cone. However, in both field models the small values of t in the southern hemisphere 

u \\ 

are found near \jjj = 50°, thus placing the non-Io C beaming pattern on the -#i side of 
the cone. No more can be said about this in the present context, and further discussion 
will be deferred to the companion paper. 

At 16 MHz, both the observed source morphology and plots of t versus Aj|j are 
I qualitatively the same as at 20 MHz. 

Before describing the beaming model at 10 MHz, we should review the following 

A. 

observations, ThC> observed source locations have shifted noticeably at 10 MHz. This 
shift has been investigated in some detail by Thieman and Smith (1978), who' suggest that 
between 15 and 10 MHz, Sources A and B have gradually moved together, merging by 

10 MHz into a single source at 210° CML. Between 22 and 10 MHz the position of 

• ■; 

Source C increases linearly, reaching a CML of about 350° at 10 MHz. The Fourth Source, 
D, near 60° CML becomes apparent below 18 MHz, though lo-independent radiation from 

fl 

that source has not been observed. 

The theoretical situation at 1 0 MHz is also more complicated than at higher fre- 

j; O 

quencies. To fit the model to the observation that there is a single source at Am as 200°, 
some change in 'h is necessary. The minimum value of 'h that still results in emission into 

the ecliptic from L ^ 6 is 4* as 40°. Emission from the northern hemisphere from L — 

-.1 

5—6 is then concentrated around Aju as 200° (Figure 6a). [Note that at 10 MHz, the O4 
and JPL models are essentially identical in all respects in both hemispheres.] With 'P as 
40°, however, the emission cone in the southern hemisphere barely intersects the 
ecliptic at all at L as 6. It is possible that very little emission from a 40° cone could 
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be beamed into the ecliptic. But if emission from the south were beamed into a cone 
with 'I' = 60° (Figure 6b), then the preferred directions (small i) would be Aju = 110° 
and 360°, close to positions of lo-D and non-lo C» respectively. One implication of this 
scenario at 1 0 MHz is that whatever the emission mechanism, more intense radiation should 
be produced as decreases from 79° to 40°, at least at low frequencies. Or, alternatively, 
radiation emitted locally at relatively large 'h and low frequency is focused by refraction 
into a smaller emission cone as the radiation leaves the Jovian environment. (At high fre- 
quencies — 20 MHz and above, Figures 4 and 5 — emission from L = 6 cannot intercept 
the ecliptic from either hemisphere if 4' is as small as 40”.) 

■ ’ / 

If one takes L S 3 as the appropriate source location of the Io-indeptfndent DAM, 

T;J) 

then 'I' s 20° would be necessary to produce a minimum in 4 at Am between 170°— 230° 
from the north, but 'F = 60° would still suffice in the south. Furthermore, radiation 
from L = 2 is beamed almost directly into the ecliptic, and should show little, if spy, de- 
peudence on CML. In fact, although the occurrence probability of the Io-independent > --==^ 
DAM (peak flux < 4 x 10^ Jy; Deschd976) is modulated in CML, there remains nearly 
a 50% occurrence probability of detecting radiation at any CML (v. Dulk and Clark 1966 
and Thieman 1977). This illustrates a systematic trend that is supported by these geo- 
metrical constructions; viz., that while the minima in t as a function of Am indicate re- 

ft ■ 

gions of preferred beam ing- -in to the ecliptic, as t decreases with decreasing frequency for a 
given L-value, an increasing fraction of the emission should become independent of CML. 

d) The Declination Effect 

It has been known for some time that the longitude of the center of Source A is 
closely correlated with Dg (v. Carr and Gulkis 1969, and Carr and Desch 1976, for reviews). 
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This sinusoidal drift of Source A is thought to occur oniy in the non-Io controlled source 

' t v ^ ! 

(Register and Smith 1969, Goertz 1971). The variation with Dg is not a simple shift of 
the center of the occurrence probability pattern, but is rather more intricate. Source A 
emission usually begins at the same CML (200°) independent of Dg. However, as Dg 
becomes more positive, the high longitude boundary of Source A moves from 280° CML 


W 


to higher longitudes, reaching 320° CML when Dg - +3.3°, a change of some 40°.' 

The centroid of the source also drifts, but by a lesser amount (s 15°) (Bozyan and Doug- 

1 

las 1976). Although Donivan and Carr (vl969) have reported finding a weaker shift of the 
Source B centroid that was in phase with the Source A drift, Thieman (1977) found, only 


a hint of an in phase drift for non-Io B, which may not have been statistically significant. 

The geometrical beaming model we have proposed for the locations of the various 

\\ 

non-lb sources can also explain this variation in Sources A and B with Dg. One of the 
effects of varying Dg is intuitively obvious, though unimportant as we shall \see. As Dg 
increases, the two sides of the emission cone that intersect the ecliptic will tend to move 
apart; i.e„ 2jt will tend to increase, producing an anti-phase variation in the centroids of 
Sources A and 53, contrary to observation. That this is actually a very small effect, an 
order of magnitude smaller than the observed 40° expansion in Source A, can be seen by 
noting that the centroids of non-Io A and B at, say, 20 MHz, are separated by some 140°, 
while the maximum allowed separation, assuming a 79° emission cone, is 158°. The 
change in 2 (i is actually much less than this maximum of 20°, for Dg changes by at most 
6.6° and so the maximum cannot be reached. With i = 45° (cf. Figures 5a and 6a) and 
using equation (3), the variation in 2ji is found to be of order 5°. The variation is even 
smaller if the radiation originates at lower L-values. 


There is, however, another explanation. For the sake of this exposition, assume 
that non-lo DAM originates at L = 6, tiiat the JPL model gives the appropriate variation 
of i with CML at 20 MHz, and 'F = 79°. [This argument cun be modified to utilize the 
0 4 model and/or a source location inside L = 6.1 When D E = -3.3°, the extent of non-lo 
A ts some 80°. in Figure 7, part of the L = 6 curve from Figure 5a is redrawn. A line, 
denoted D E = -3.3°, at < ss 44° intersects the Source A region at points (a) and (b) sepa- 
rated by 80°. Assume that above the lino, radiation is received in the ecliptic, but radia- 
tion emitted below the line with larger t is not. When Dg is increased to +3.3°. a second 
line can be drawn at i =* 50°, which intersects the Source A region at points (efi'md (d). 
The high longitude boundary of Source A has moved some 50°. while the beginning of 
Source A lias moved at most by 25° (probably less because of potential interference from 


Source B). 


Furthermore, the longitudes of the Source A and B centroids will increase in phase. 


This argument is only meant to be suggestive. Our present knowledge of the near- 
surface Jovian magnetic field is too imprecise for more detailed modeling. Furthermore, 
refraction and other propagation effects arc likely to play an important role in any com- 
plete explanation of this phenomenon. However, several general features are germane. 
First, the fact that the high longitude boundaries of Sources A and B may move in phase 
by as much as 40° does not necessarily contradict predictions based on Dulk’s (1967) 
conical sheet model, since the motion proposed here is not caused by a tilting of the 

cone axis. Second, there is amplib reason to expect that inhomogenoities in the neur sur- 

V 

face Jovian fields will produce the very asymmetrical drifts of the low and high longitude 
sides of Source A that are observed, and that we have described at: least qualitatively in 
terms of this beaming model. Two caveats should be noted. One is that we cannot esti- 

• *■ -y 

mate how much the centroids of non-Jo A and B will drift, except to say that for non-lo 
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A it should be nibre than the 5° that results from the change in 2 /j discussed above, >r 
non-lo B, thedrii't due to changes in 2fi will tend to cancel the change in the centroid 
produced by the asymmetrical expansion of B to higher longitudes. The expected result 
is that the drift of the centroid of B will be less than that of A, as observed. The second 
is that, although this construction implies' iiiat the low longitude side of A (and B) should 
drift much less than, and in a direction opposite to, the high longitude side, there is no 
observational evidence to suggest that the low longitude side of Source A (or B) drifts at all. 
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111. CONCLUSIONS 

The theoretical model we have presented is, we believe, useful for explaining some 

Q ! - 

of the significant features of the lo-independent Jovian decameter radiation. 

* It is based on four assumptions, which hopefully will be tested experimentally on 
future space probes to Jupiter. The first of these is that the Io-independent DAM is pro- 
duced by electrons stablyl^pped in Jupiter’s inner magnetosphere (L ^ 6); second, that 
the radiation is excited near the local gyrofrequency in the R-X mode; third, that the 
Q radiation is beamed into a conical surface; and last, that the beaming of radiation into 
the ecliptic depends sensitively on the magnitude of the till of the emission cone toward 


the Jovian equatorial plane. 

5a==s rh$ first two assumptions imply that the Io-independent sources radiate below 
29 MHz, for only precipitating particles can reach higher gyrofrequencies. Thus we do 
not expect Io-independent radiation to be seen much above 29 MH/.. 


Above 10 MHz, this model suggests that decameter radiation is beamed into a 

conical surface with a half-angle, that is between 60°— 80°. The lower value obfijins 

*> 

if radiation originates at small L-values (L 3^3), while emission near L 2= 6 requires 4' s 

'N . 

80°. At 10 MHz, it was necessary to choose 4' Independently in the twq hemispheres — 

40° in the north, and 60° in the south. 

' \! ■ . » ■ 

The model,. re constructed is then consistent with several features off PAM morphol- 
ogy, The association of lew values of t with high occurrence probabilities leads to the 

* ft ■ 

conclusion that at 27 MHz, radiation is preferentially beamed into the ecliptic from only 
the northern hemisphere, and only at the central meridian longitudes of the non-lo A and 
B sources. Below 27 MHz, two additional regions of small i become apparent, originating 
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in the southern hemisphere. One of these lies near the location of non*Io C. But, as 
mentioned above, there are difficulties with this interpretation. 

//At J O MHz, with suitable variation in ' 4 ', the beaming pattern shifts so that there 
is a single source from the northern hemisphere, near the observed location (200° CML), 
while two sources remain in the south, one of which may correspond to non*Io C. Be- 
cause of the generally low values of t at all longitudes, the model is consistent with the 
observed tendency at low frequencies for there to be considerable DAM noise independent 
of CML. ' 


An additional success of the model is the suggested explanation of the large drift 
of the high longitude boundary of non-Io a with changes in Dg. 

\7 

The model also suggests the possible existence of DAM features that have not yet 

S \| 

been observed. These include a second southern hemisphere (left-circularly polarized) 
source near 150° CML (f > 10 MHz), as well as Jo-independent emission from 120° 


CML (Sourer, D) at 10 MHz. Because of the sensitrrity of the beaming pattern to small 
changes in t, it is also expected that when Jupiter is observed by space probes at high '' 


northern and/or southern latitudes, two trends will appear. One is an increased occurrence 


I 
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probability at high frequencies (f =5- 27 MHz.) that is still strongly modulated with longitude. 
This is expected, because although one will be well inside the emission cone as it inter- 
cepts the receiver, the restricted azimuthal region within which, say, f^ > 27 MHz will 

\> 

continue to produce a strong modulation with longitude. The second effect, most notice- 
able as the frequency decreases below 16 MHz, would be a decline in the modulation of 
P with CML as the emitting region gradually expands to encompass nearly all longitudes. 


One limitation of this geometrical model is the absence of an explanation for the 
significant differences in occurrence probability among the various lo-independent sources. 
For example, non-lo A generally has a higher occurrence probability than non-Io B (Carr 
and Desch 1976). The explanation for these differences may well depend on exactly how 
the inhomogeneities in the local magnetic fields affect both the amplification and propa- 
gation of the radiation in the source region. In general, these strong gradients in the near 
surface magnetic fields, where the high order multipole terms are very important, will not 
be symmetric about the emission cone. This is suggested by even cursory examination of 
Figures 3 and 2, and can be confirmed by more detailed calculations using either magne- 
tic Held model. Similarly, the effects of these gradients in the magnetic field cn an elec- 
tromagnetic wa'vc propagating through this underdense plasma (o> pe « |S2 e |) should ex- 
ceed the effects of any gradients in the density at these altitudes. Thus, refraction will 
tend to bend the wave number vector, k, toward the direction of -3£nB/3r (Stix 1962, 

/v 

Smith 1973). 

- 

The maximum gradient in B at these low altitudes lies close to the direction of B 
itself. Consequently, a ray initially launched at ’if 80° will be bent toward the axis of 
tile emission cone. Unless t is small (just small is difficult to estimate without doing 
a complete ray-tracing calculation, which is beyond the scope of this paper), refraction 
toward the cone axis will have the effect of diminishing the amount of radiation that is 
beamed into the ecliptic. Because the gradients in B are asymmetric about the emission 

K /w 

cone axis, the amount of radiation beamcd\ into the ecliptic should differ on the two sides 

•' j! jf 

of the cone, at least partially accounting for the differing occurrence probabilities of the 
various sources. ; V 
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■(} In the accompanying paper, the features of this geometrical model are combined 

with a specific instability mechanism in an attempt to reduce the number of assumptions 

needed to describe the Io-independent DAM. That analysis retains the assumption of 

excitation of DAM by trapped electrons. It is then found that ihe value Of >]> and the 
■ // - . ■ 
cone thickness, A 'I', can be determined as a function of frequency, and that the L-shel! 

of the source region can be localized within some fairly well defined Limits. The theory 

is also able to account for the excellent correlation in the northern hemisphere between 

-_.fi 

small. values of t and the source locations. A solution to the difficulties encountered in 

X 

explaining the geometry of the southern hemisphere source is also suggested, in which 

'* ~ 'x • ' 

both the Io and non-lo C sources are emitted into the +fi side of the emission cone. 
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mirror field exceeds 27, 20, 16 or 10 MHz (shaded areas and dashed line). The contours are projected onto the alti- 
tude of the cloud tops. Also drawn are isointensity contours in 1G increments of Jupiter’s surface field. Footprints on 
the surface are plotted for L = 6 (“+”), L = 5 (“X”), L = 4 (“*”), L = 3 (“0’i') and L = 2 (“#”) field lines. The azimuthal 
spacing of the field lines in the equatorial plane is 10°. Conjugate points at aiiiy longitude can be found by noting that in 
the two hemispheres, near Xjjj = 0°, footprints denoted “S” are conjugate. (CI4 model). „ 


SOUTH 



Figure 2. Similar to Figure 1, but using the JPL model 



Figure 3. A plot of i versus Aju at 27 MHz, using the O 4 model with = 79°. Gaps in the 

ciirves.j'pr, L = 3 and 5 result from localized regions where fjyj is less than 27 MHz. Curves on 
the right are solutions of Ajjj - Xj jj + p while those to the left are obtained from Ajjj = X u j 
fi. The curves represent the central meridian longitude directions of the intersection lines of 
the emission cone with the ecliptic. Note that there is no longitude for L = 2 at which 
f M ^ 27 MHz (cf. Figure 1). Preferential beaming into the ecliptic t- small values of t — 
occurs at the positions of the Io independent Early and Main Sources (3 and A, respectively). 
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0 4 MpOEL (SOUTH) 



Figure 6. A plot of i versus Ajjj at 10 MHz, using the 0 4 model. Use of the JPL model produces identical curves. In the north 
(Figure 6a) 'P = 40°, in the south (Figure 6b) M> = 60°. See the text for a discussion of the source locations. 





JPL MODEL (NORTH) = 79 


20 MHz 


Figure 7. A possible explanation of the Declination Effect is shown. As Dg changes from 
-3,3° to +3.3° , the width of non-Io A can increase from AAjjj = (b) - (a) to AAjh - (d) - 
(c). (JPL model on L = 6 , northern hemisphere). 







FIGURE CAPTIONS 


Figure L Two dimensional contour plots showing regions within which, f M > the maximum 
electron gyrofrequency at the mirror field B^, exceeds 27, 20, 16 or 10 MHz (shaded 
areas and dashed line). The contours are projected onto the altitude of the cloud 
tops. Also drawn are isointensity contours in 1G increments of Jupiter’s surface field. 
Footprints on the surface are plotted for L = 6 (“+”), L = 5 (“X”), L - 4 (“*”), 


L = 3 (“0”) and L = 2 (“#”) field lines. The azimuthal spacing of the field lines in 
the equatorial plane is 10°. Conjugate points at any longitude can be found by noting 
that in the two hemispheres, near Xjjj = 0°, Footprints denoted “S’* are conjugate. 

4 

(O4 mocjiel). ''' 

Figure 2.i Similar to Figure 1, but using the JPL model. 


J 

Figure 3t>: A plot of t versus Ajjj at 27 MHz, using the 0 4 model with 'P - 79°. Gaps in the 
curves for L = 3 and 5 result from/iocalized regions whsre f M is less than 27 MHz. 

f 

Curves on the right are solutions of A lt j = Xjjj + fi, while those to the left are obtained 
from Am - Xjjj - p. The curves represent the central meridian longitude directions of 
the intersection lines of the emission cone with the ecliptic. Note that there is no 
longitude for L = 2 at which fj^ > 27 MHz (ef. Figure 1). Preferential beaming into 
r the ecliptic — small values of t — occurs at the positions of the lo independent Early 
and Main Sources (B and A, respectively). 

Figure 4. Similar to Figure 3, but for 20 MHz. For clarity, curves are not plotted for L = 4 
or 5 ; and, northern and southern hemispheres are plotted separately (Figures 4a and 
4b, respectively). In the north, minima in 1 (peaks on the curves), occur near the 
CML of the non-Io A and B sources. Note that at L - 6, the preferred emission region 
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for non-Io A extends to Third Source (C) longitudes. See the text for a discussion 
of the southern hemisphere sources. 

0 

Figure 5. Similar to Figure 4, but using the JPL model. Note that for L = 2, the allowed 

range of Ajjj is much less than it was using the 0 4 model (cf. Figure 3). In addition, 
at L = 6 in the north (Figure 5a), small values of i do not extend to non-Io C longi- 
° tudes. See the text for a discussion of the southern hemisphere source locations 
(Figure 5b). 

Figure 6. A plot of t versus Am at 1 0 MHz, using the 0 4 model, Use of the JPL model pro- 


duces identical curves. In the north (Figure 6a) 'P - 40°, in the south (Figure 6b) 
^ = 60°. See the text for a discussion of the source locations. 


Figure 7. A possible explanation of the Declination Effect is shown. As Dg changes from 



-3.3° to +3.3°, the width of non-Io A can increase from AAjj[ = (b) - (a) to A Am = 
(d) - (c). (JPL model on L - 6, northern hemisphere). 
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